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ABSTRACT 

We derive an integral condition for core-collapse supernova (CCSN) explosions and use it to construct 
a new diagnostic of explodability. The fundamental challenge in CCSN theory is to explain how 
a stalled accretion shock revives to explode a star. In this manuscript, we assume that the shock 
revival is initiated by the delayed-neutrino mechanism and derive an integral condition for spherically 
symmetric shock expansion, v s > 0. One of the most useful one-dimensional explosion conditions is 
the neutrino luminosity and mass-accretion rate (L„- J\4) critical curve. Below this curve, steady-state 
stalled solutions exist, but above this curve, there are no stalled solutions. Burrows & Goshy (1993) 
suggested that the solutions above this curve are dynamic and explosive. In this manuscript, we 
take one step closer to proving this supposition; we show that all steady solutions above this curve 
have Vs > 0. Assuming that these steady v s > 0 solutions correspond to explosion, we present a 
new dimensionless integral condition for explosion, 'F > 0. T roughly describes the balance between 
pressure and gravity, and we show that this parameter is equivalent to the r condition used to infer the 
L v -M critical curve. The illuminating difference is that there is a direct relationship between 'P and 
v s - Below the critical curve, *P may be negative, positive, and zero, which corresponds to receding, 
expanding, and stalled-shock solutions. At the critical curve, the minimum 'll solution is zero; above 
the critical curve, ’F m in > 0, and all steady solutions have v s > 0. Using one-dimensional simulations, 
we confirm our primary assumptions and verify that ’Pmm > 0 is a reliable and accurate explosion 
diagnostic. 

Keywords: supernovae: general hydrodynamics — methods:analytical methods: numerical 
shock waves 


1. INTRODUCTION 

The question of how typical massive stars explode 
as core-collapse supernovae (CCSNe) has plagued theo¬ 
rists for decades (Colgate & White 1966). The collapse 
and bounce of the Fe core of a massive star launches a 
strong shock wave; however, this prompt shock quickly 
stalls as a result of electron capture, nuclear disso¬ 
ciation, and neutrino emission (Hillebrandt & Mueller 
1981; Mazurek ct al. 1982; Mazurek 1982). Thus the 
explosive shock is momentarily aborted, producing a 
stalled accretion shock. If the shock remains stalled, 
the explosion will fail and the proto-neutron star will 
continue to accrete, eventually collapsing to a black 
hole (Fischer et al. 2009; O’Connor & Ott 2011). We 
know, however, that stars explode and many times 
leave neutron stars (Li et al. 2010; Horiuchi et al. 2011; 
Fryer et al. 2012). Therefore, the fundamental question 
in core-collapse theory is how a CCSN transitions from 
a stalled accretion shock phase into a phase of runaway 
shock expansion that explodes the star. In this paper, 
we derive an integral condition in the limiting case of 


spherical symmetry that divides stalled-shock solutions 
from explosive solutions. 

The prevailing view is that neutrinos help to reinvig¬ 
orate the stalled shock, leading to runaway expansion 
(Bethe & Wilson 1985; Janka 2012; Burrows 2013). 
A large neutrino flux cools the proto-neutron star 
and surrounding regions; about half of that luminosity 
comes from neutrinos that diffuse out of the natal proto¬ 
neutron star, and the rest is emitted directly as accretion 
luminosity by the accreting stellar material. A fraction 
of these neutrinos recapture just below the shock, 
depositing enough energy in the post-shock mantle to 
reinvigorate the stalled accretion shock and initiate 
explosion. However, when this picture is simulated in 
numerical models with significant (although perhaps not 
sufficient) detail, the final and most important element 
explosion — remains elusive or, at best, inconsistent 
with observations (Ott et al. 2008; Muller et al. 
2012b; Hanke et al. 2013; Takiwaki et al. 2014; 
Lentz et al. 2015; Dolence et al. 2015; Melson et al. 
2015; Brucnn ct al. 2016; Roberts et al. 2016). There- 
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fore, a major challenge for simulators has been to 
understand the various physical effects that influence 
the unsatisfactory outcomes, all in the context of 
enormously complicated and expensive simulations. 
This process of interpretation has inevitably relied 
on some easily understood models that, one hopes, 
encapsulate the important physics of the problem. 

In this context, perhaps the most impactful model 
was proposed by Burrows & Goshy (1993). They pro¬ 
posed that the essence of the core-collapse problem 
is captured by considering a simple boundary value 
problem. The inner boundary is set at the proto¬ 
neutron star “surface,” defined by a neutrino optical 
depth of 2/3, where a temperature and therefore lu¬ 
minosity is specified. The outer boundary is at the 
stalled shock, where the post-shock solutions match 
the nearly free-falling stellar material upstream via the 
jump conditions. Burrows & Goshy (1993) proceeded 
to show that hydrodynamic solutions with a station¬ 
ary shock only exist below a critical luminosity and 
accretion-rate curve. Above this curve, no stalled- 
shock solutions exist. They speculated that this critical 
curve, separating stationary from non-stationary solu¬ 
tions, also represents a critical curve for explosion; they 
suggested (but did not prove) that non-stationary so¬ 
lutions are explosive. Subsequent work with parame¬ 
terized one-, two-, and three-dimensional core-collapse 
models have confirmed qualitatively that such a criti¬ 
cal curve exists (Murphy & Burrows 2008; Hanke et al. 
2012; Dolence et al. 2013; Couch 2013), but its precise 
location in parameter space depends on other details of 
the problem, including, e.g., the structure of the pro¬ 
genitor, which hampers the use of the critical curve in 
practice (Suwa et al. 2014; Dolence et al. 2015). 

Nonetheless, the idea of criticality has framed 
much of the discussion surrounding the simula¬ 
tion results and has motivated the introduction of 
heuristic and approximate measures of “nearness to 
explosion.” For example, many have suggested 
that the ratio of the advection timescale to heat¬ 
ing timescale captures an important aspect of the 
problem, with values >1 conducive to explosion 
(Janka & Keil 1998; Thompson 2000; Thompson et al. 
2005; Buras et al. 2006a; Murphy & Burrows 2008). 
Similarly, Pejcha & Thompson (2012) argue that their 
“antesonic” condition represents a critical condition for 
explosion. Both conditions, however, suffer the same 
afflictions: they lack precise critical values, and they 
both run away only after explosion commences. These 

1 Time to advect through the net heating region 

2 Time to significantly change the thermal energy in the net 
heating region 


shortcomings have led to the practice of measuring these 
parameters as a function of simulation time and de¬ 
ciding on “critical values” ex post facto (Muller et al. 
2012b; Dolence et al. 2013; Couch 2013). Clearly, this 
is unsatisfactory and leads to little more insight than 
simply identifying explosions with runaway shock radii. 
In principle, these critical values may even be reached 
precisely because of explosion; the relationship may be 
symptomatic rather than causal. 

Investigating criticality has proven to be a useful 
quantitative measure of explodability. For some time, it 
was very clear that the delayed-neutrino mechanism fails 
in one-dimensional, spherically symmetric simulations 
(Liebendorfer et al. 2001b, a, 2005; Rampp & Janka 
2002; Buras et al. 2003; Thompson et al. 2003); it was 
also becoming apparent that multi-dimensional simu¬ 
lations showed promise where one-dimensional simula¬ 
tions failed (Herant et al. 1994; Janka & Muller 1995, 
1996; Burrows et al. 1995, 2007; Melson et al. 2015; 
Bruenn et al. 2016; Roberts et al. 2016). The critical 
curve offers a quantitative measure of how much the 
multidimensional instabilities aide the neutrino mecha¬ 
nism toward explosion. Murphy & Burrows (2008) in¬ 
vestigated in one and two dimensions whether a crit¬ 
ical curve is even sensible in time-dependent simula¬ 
tions. They not only empirically found that a criti¬ 
cal curve is viable, but that it is about 30% lower in 
two-dimensional simulations. Subsequent investigations 
indicate that the critical curve in three-dimensional is 
similar to two-dimensional simulations (Couch 2012; 
Hanke et al. 2013). There are certainly slight differences 
between two and three dimensions, but these differences 
are minor in comparison to the major shift that mul¬ 
tidimensional instabilities enable in going from one di¬ 
mension to multidimensionality. Further investigations 
suggest that turbulence plays a major role in reducing 
the critical condition for explosion (Murphy fc Meakin 
2011; Murphy et al. 2013). While multidimensionality 
and turbulence are important considerations in the ex¬ 
plosion mechanism, we do not address turbulence in this 
paper; we leave that for a subsequent paper (Q. Mabanta 
et al., in preparation). In this manuscript, we aim to fur¬ 
ther understand the foundational aspects of criticality, 
and later, we hope to expand these to include turbu¬ 
lence. 

In this work, we revisit the idea of criticality, gener¬ 
alizing the critical-curve concept to a critical hypersur¬ 
face that depends on all of the relevant parameters. We 
show that for a fixed set of parameters, there is a family 
of possible solutions. Depending upon the parameters, 
each family falls into one of two categories. In one cate¬ 
gory, the family consists of solutions with negative, zero, 
or positive shock velocity. For such a family, the zero- 
shock-velocity solution is the quasi-stationary solution. 
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In the second category, all possible solutions have pos¬ 
itive shock velocity. These two categories are divided 
by a critical hypersurface, which may be expressed as a 
single dimensionless parameter. 

Our initial motivation in revisiting criticality was to 
ask why does the neutrino-luminosity and accretion-rate 
curve of Burrows & Goshy (1993) corresponds to explo¬ 
sion in simulations. Under certain restrictive but empir¬ 
ically reasonable assumptions, we derive that the only 
possible solutions above the curve correspond to positive 
shock velocity. Not only does this derivation suggest a 
reason for explosion, but it also suggests a more general 
critical condition. 

This new condition for explosion, the critical hyper¬ 
surface and associated dimensionless parameter, proves 
to be a useful diagnostic for core-collapse simulations. 
For one, we empirically show that the hypersurface cor¬ 
responds to the transition to explosion in parameterized 
one-dimensional CCSN models. Further, the associated 
dimensionless parameter reliably and quantitatively in¬ 
dicates when explosions commence. Since we derive the 
parameter directly from the equations of hydrodynam¬ 
ics, its usefulness is not limited by the ad hoc calibra¬ 
tions of other popular measures. Finally, when com¬ 
bined with semi-analytic models, we show that our sin¬ 
gle parameter yields an accurate and reliable measure of 
nearness-to-explosion. 

For your convenience, the structure of this manuscript 
is as follows. In section 2, we review the boundary 
value problem that describes the core-collapse problem 
(Burrows & Goshy 1993) and identify the important pa¬ 
rameters of the problem. With this framework, and with 
the proposition that v s > 0 corresponds to explosion, we 
derive the integral condition for explosion in section 3. 
Then in section 4, we validate the integral conditions 
for explosion with one-dimensional parameterized sim¬ 
ulations. In section 5, we use the integral condition 
to investigate the family of steady-state solutions and 
propose an explosion diagnostic. Then in section 6, we 
compare the reliability of the integral condition explo¬ 
sion diagnostic with other popular explosion measures, 
finding that the integral condition outperforms them all. 
Finally, for a brief summary of the conclusions, caution¬ 
ary notes, and future prospects, see section 7. 

2. QUASI-STEADY SOLUTIONS: A BOUNDARY 
VALUE PROBLEM 

The fundamental question of core-collapse theory is 
how does the stalled accretion shock phase transitions 
into a dynamic explosion. In other words, what are the 
conditions for which the shock velocity, v s , is persis¬ 
tently greater than zero? While the shock is stationary, 
the steady-state assumption is quite good. Under this 
assumption, the entire region below the shock may be 


treated as a boundary value problem (Burrows & Goshy 
1993). The upper boundary is set by the properties of 
the nearly free-falling stellar material and the jump con¬ 
ditions at the shock, while the lower boundary is the sur¬ 
face of the neutron star. Later, we derive a condition in 
which there are no more v s = 0 solutions, and the only 
solutions left are those in which v s > 0. To understand 
when the steady-state solutions are no longer viable, we 
must first understand the stalled accretion shock solu¬ 
tion and the important parameters of the problem. 

To begin, the governing conservation equations are 

^+V-(pv) = 0, (1) 


and 


d(pv) 

dt 


+ V • (pvv) + VP 


pS7(t>, 


d(pE) 

dt 


V- 


pv 



pv- V0 + pg, 


( 2 ) 

( 3 ) 


where p is the mass density, v is the velocity, P is the 
pressure, <fi is the gravitational potential, E = e+v 2 /2 is 
the internal plus kinetic specific energy, and h = e + P/p 
is the specific enthalpy. In this paper, we approximate 
the heating and cooling, q as 


_ L u k 
q 4-7T r 2 


- C 



( 4 ) 


The first term is a model in which we treat the neutrino 
heating as if all neutrinos were emitted from the proto¬ 
neutron star “surface” with a luminosity of L v ( anka 
2001; Murphy & Burrows 2008; Murphy et al. 2013). n 
is the opacity for absorbing neutrinos in the region above 
the “surface,” which is proportional to T 2 . In this pa¬ 
per, the surface of the proto-neutron star is given by 
r = J^° up dr = 2/3, which implicitly defines the proto¬ 
neutron star radius (R v ) as the neutrinospheric radius. 
In practice, we find that this corresponds to a density of 
P 2/3 ~ 7 x 10 10 g cc -1 . The second term is the neutrino 
cooling due to thermal weak interactions (Janka 2001). 

The neutrino interactions that we employ are quite 
simple, and while they are crude, they retain some key 
elements allowing us to more easily assess the impor¬ 
tance and role of neutrino heating and cooling. In the 
simple model, there are three parameters to describe 
the neutrino heating, the core neutrino luminosity, the 
temperature of the neutrinos, and the neutrino-sphere 
radius. A better description of the neutrino luminos¬ 
ity would include an accretion luminosity, which would 
account for the added neutrino luminosity provided by 
the semi-transparent cooling region above the neutrino 
sphere. While it is not difficult to include such an accre¬ 
tion luminosity (Pejcha & Thompson 2012), we decided 
to employ the more traditional “light bulb” description 
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to reduce the confounding prescriptions in our initial 
investigation. Later, we will include the accretion lumi¬ 
nosity and investigate the differences and whether they 
are merely quantitative or are more substantially quali¬ 
tative. Even though the neutrino spectrum is not quite 
Planckian, we treat it as such to reduce the neutrino en¬ 
ergy spectrum to one parameter, the neutrino tempera¬ 
ture. The final parameter, the neutrino sphere radius, is 
probably even more approximate. Even if one would be 
able to describe the neutrino spectrum with one tem¬ 
perature, there is a distribution of neutrino energies, 
and since the cross section depends upon the square of 
the neutrino energy, the neutrino sphere for each energy 
would occur at a different radius. Nonetheless, we seek 
a parameterization of the size of the proto-neutron star 
and choosing one neutrino sphere radius which emits 
the core luminosity at a core temperature probably pre¬ 
serves the qualitative behavior. Given these differences, 
we present the ideas and results of this paper as qual¬ 
itative arguments. We suspect that subsequent more 
realistic models for neutrino interactions will preserve 
the qualitative nature of our conclusions. 

In a quasi-steady-state, the time-derivative terms are 
small and the problem is well-approximated by setting 
these terms to zero. We specifically refer to the quasi¬ 
steady-state because we are assuming that the post¬ 
shock profile is steady, but the shock velocity is nonzero. 
Often when one assumes steady-state, one also assumes 
that the shock velocity is zero. However, it is possible to 
have a quasi-steady nonlinear solution behind the shock 
and a nonzero shock velocity. The Noh test problem is 
an example of one such nonlinear solution. In this test 
problem, a supersonic flow is incident on a wall. A shock 
forms, and the post-shock flow maintains the same den¬ 
sity profile; i.e. it is steady. However, the shock has a 
nonzero lab frame velocity. Similarly, in finding quasi¬ 
steady solutions, we assume that the post-shock flow is 
steady (i.e. the time-derivative terms are zero), but we 
explicitly allow for a nonzero shock velocity. 

The solution to the resulting quasi-steady boundary 
value problem describes the conditions of the flow be¬ 
tween the proto-neutron star surface and the bounding 
shock in terms of the important parameters of the prob¬ 
lem: L„, R v , Mns (proto-neutron star mass), and 

M (accretion rate). Even though we have assumed that 
the post-shock flow is steady, the shock velocity may 
or may not be zero. We ask the fundamental question 
what it takes for all quasi-steady-state solutions to have 
v a > 0. Later, we will show that for certain values of 
these important parameters all of the quasi-steady state 
solutions have v s > 0. See Figure 1 for a schematic 
of the boundary value problem and the most important 
parameters of the problem. 

The five parameters that we highlight are a natural 


parameterization for the core-collapse problem and of¬ 
fer a way to parameterize some of the most uncertain 
aspects of the core-collapse problem. For example, two 
uncertain aspects of the core-collapse problem are the 
structure of the progenitor and the dense nuclear equa¬ 
tion of state (EOS). 

The uncertainties of the dense nuclear EOS for neu¬ 
tron stars are often parameterized in terms of a mass- 
radius relationship (Lattimer & Prakash 2016). Each 
EOS predicts a specific mass-radius curve. There¬ 
fore, investigating how the conditions for explosion de¬ 
pend upon the proto-neutron star mass, Mns, and the 
neutrino-sphere radius (a proxy for the neutron star ra¬ 
dius), R v , provides a means to parameterize how the 
condition for explosion depends upon the uncertainties 
in the EOS. 

The three other parameters help to parameterize the 
uncertainties in the progenitor structure as well as neu¬ 
trino diffusion in the core. For example, the time evolu¬ 
tion of both Mns and M depends upon the progenitor 
structure. L„ depends upon a mix of the total thermal 
energy available in the neutron structure, neutrino dif¬ 
fusion, and the accretion rate. Hence, L v is not entirely 
independent of the other parameters. In a subsequent 
paper, we will explore the consequences of these extra 
constraints. For now, however, we use L u to discuss 
our new explosion condition in the context of previous 
conditions. The neutrino temperature is also not en¬ 
tirely independent (e.g. L u ss AnR^aT^), but again 
we use it to make a connection to past literature. In 
a future paper, in which we discuss analytic solutions, 
we will propose an alternative formulation for these pa¬ 
rameters. Until then, these particular parameters are 
a natural parameterization of the core-collapse problem 
that also incorporate the uncertainties in the progenitor 
and neutron star physics. 

3. DERIVING AN INTEGRAL CONDITION FOR 

Us >0 

We reduce the core-collapse problem to a set of in¬ 
tegral conditions and show that this leads to a critical 
condition for explosions. Before we derive the integral 
conditions, let us describe what motivated us to con¬ 
sider the integral equations at all. The governing equa¬ 
tions (eqs. 1-3) are commonly presented in differential 
form, but one could just as easily present them in in¬ 
tegral form. The two forms represent exactly the same 
information. The choice is determined by the ease and 
method for finding solutions. For example, it is typically 
easier to work with the differential form when one needs 
to find numerical solutions, but if one seeks an analytic 
solution, the integral equations can be easier to use. For 
example, in finding the motion of a body in a potential 
field, one may solve the equations of motion, or one may 
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What are conditions for 
v s = 0 to become v s > 0? 


Figure 1. After the bounce shock stalls, the post-shock flow settles into a quasi-steady configuration. The solutions satisfy a 
boundary value problem between the neutron star “surface” and the stalled shock. The important parameters in determining 
the steady-state solutions are the neutrino luminosity being emitted from the core (L v ), the temperature of the neutrinos (T„), 
the mass-accretion rate onto the shock and NS (Af), the neutrino sphere or NS radius {R v ), and the NS mass (Mns)- The 
fundamental challenge in core-collapse theory is to understand how v 3 = 0 transitions to v s > 0. In section j, we derive such a 
condition. 


use the integral condition, conservation of energy. 

We argue that the integral equations provide an eas¬ 
ier route to deriving a unified condition for explosions. 
Because the shock is a crucial component of the core¬ 
collapse problem, when we derive the explosion condi¬ 
tion we embark on a route that is similar to deriving the 
Rankine-Hugoniot shock jump conditions. The route to 
deriving these jump conditions involves using the inte¬ 
gral equations in the context of a moving boundary, the 
shock. 

3.1. The generic integral condition 

As in deriving the Rankine-Hugoniot conditions, we 
start with the conservation equations and the Reynolds 
transport theorem (or Leibniz-Reynolds transport the¬ 
orem). Consider a generic conserved quantity, q , which 
could represent a mass density, momentum density, en¬ 
ergy density, etc. Generically, one may define the to¬ 
tal conserved quantity in region D(t) as J D ^ qdV and 
the time-rate-of-change of this conserved quantity is 
It (fo(t) Q . The Reynolds transport theorem de¬ 
composes this time-rate-of-change into an Eulerian com¬ 
ponent and a component that accounts for the motion 


of the boundaries of region D(t): 

-tA i < l dV \ = [ ~T: dV + / qVb-dS , ( 5 ) 
at \JD(t) J JD(t) at JdD(t) 

with dD(t) the surface of this domain, dS the surface 
element, and Vf, the velocity of the surface. To make 
use of the Reynolds transport theorem, we need an ex¬ 
pression for the Eulerian time derivative, dq/dt. The 
conservation equation provides such an expression. A 
general conservation equation has the following form: 


where f is the flux of q , and s is a source/sink of q. With 
these two equations (eqs. 5 & 6) we may now proceed 
to derive the integral conditions. 

To construct the appropriate boundary value problem, 
we must choose the appropriate boundaries. We con¬ 
sider a boundary value problem where the lower bound¬ 
ary is at r\ and the upper boundary is at r 2 . The overall 
problem is divided into two domains, with a different 
solution in each domain. The pre-shock solution is in 
the domain (r s (f),r 2 ], the post-shock solution is in the 
domain [ri,r s (f)), and the two solutions are joined by 
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the Rankine Hugoniot jump conditions at r s (t). With 
the boundaries defined, we now define the integral equa¬ 
tions. 

First, we integrate eq. (6) from n to r 2 , and since the 
boundaries are constant in time, we have 


d_ 

dt 



dq 

dt 


dV = 


(-V-f + s) dV. 


(7) 

Next, we rewrite the left-hand side of eq. (7) to explicitly 
express the shock radius. In order to do this, we must 
divide the integral into two regions, above and below 
the shock, and we use the Reynolds transport equation 
(eq. 5). This gives 


d 

dt 



= v s 4nr 2 (q s _ e - q s+e ) 



dq 

i dV+ 




( 8 ) 


where s — e and s + e represent values just below and 
above the shock. From here on, we introduce a more 
compact nomenclature for the values below and above 
the shock with and q+. 

In the core-collapse problem, the pre-shock solution is 
nearly in free-fall and can be solved analytically. There¬ 
fore, the problem reduces to finding the post-shock so¬ 
lution subject to the lower boundary, the surface of the 
NS, and the upper boundary, the shock. Now, we take 
the limit that r^ approaches r s (t). In this limit, we have 
the following simplifications: 


lim q 2 = q+ and 
r 2 


rr 2 

lim / 

r 2 ^r s (t) J rB (t) 


dq 

dt 


dV = 0, 


(9) 


which, when inserted into eq. (8), results in the simpler 
integral equation 

( 10 ) 

By plugging eq. (10) into eq. (7), we arrive at the general 
integral equation 

rsid) f) n rr a 

v s 4nr 2 (q_-q + )+J dV = 4n(f 1 rj-f+r 2 )+j s dV . 

( 11 ) 

Note that if we take the limit that rq approaches r s 
from below, this reduces to the shock jump conditions, 
but if we retain r\ at the neutron star surface, we have 
the general integral equations describing the post-shock 
structure. If we assume a steady profile of q(r) between 
the neutron star, r±, and the shock, r s , then this reduces 
to 

v s 4irr 2 s (q- - q+) = 4tt( fpr\ - f+r 2 ) + [ sdV. (12) 


Eq. (12) represents the integral condition that relates 
the shock velocity to the steady-state integrals. For the 
core-collapse problem, there is not one constraint, but 
three that come from the three equations of hydrody¬ 
namics. In principal, this could be a messy algebraic 
problem to relate v s to the integral conditions. How¬ 
ever, in the next section, we show that, in the context 
of the core-collapse problem, a simple relation emerges 
between the shock velocity and the integral conditions. 


3.2. Deriving an expression for v s in terms of the 
integral conditions 

Now, to derive the simple relation for the shock veloc¬ 
ity in terms of the integral condition. We begin by sum¬ 
marizing the traditional way of deriving the shock ve¬ 
locity in terms of the density jump and pressure jump. 
Then, we use a similar technique, but instead of the 
pressure jump, we consider the full momentum integral. 
The Rankine-Hugoniot jump conditions in mass and mo¬ 
mentum are 

p-u- = p + u+ (13) 


and 

P-+ p-u 2 _ = P++ p+u 2 + , (14) 


where once again p is the density and P is the pressure. 
This time, however, u = v — v s represents the velocity in 
the shock frame, which, pragmatically, is just the differ¬ 
ence in the lab frame velocity (v) and the shock velocity. 
As before, + denotes the pre-shock state and — denotes 
the post-shock state. Combining these two equations 
leads to the following expression for the shock velocity: 


v s = 


v s 

v+ 


-1 + 


P(P- ~ P+) 

{fi - 1 )p+v\ ' 


(15) 


where v s is the dimensionless shock velocity, scaled by 
the pre-shock speed, v + = — v+, and /3 = p~/p+ is the 
shock compression ratio. In deriving eq. (15), there are 
formally two solutions, one with + and one with — in 
front of the radical sign. One of these solutions is real 
and the other is unphysical. The solution associated 
with — corresponds to P- — P + <0 when p_ — p + >0; 
this would result in a discontinuous rarefaction wave, 
which is unphysical. The other solution with + in front 
of the radical sign corresponds to P- — P+ > 0 when 
P- — p+ > 0; this is a shock, the only physical solution. 

So far, our derivation of eq. (15) is merely a sum¬ 
mary and does not incorporate the integral conditions; 
now we use a similar procedure to derive an integral 
condition for the shock velocity. In deriving eq. (15), 


3 The jump condition in the energy equation closes the system 
by providing an expression relating the jump in density with the 
jump in pressure. 
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we began with the mass and momentum jump condi¬ 
tions across the shock; to derive the integral equivalent, 
we begin with the mass and momentum integral equa¬ 
tions for the shock velocity. Using the generic integral 
equation, eq. (12), and the mass equation, eq. (1), the 
integral equation for mass that includes the shock is 

v a r 2 8 {p- - p+) = Pivirl - p+v+r 2 ( 16 ) 

where the subscript 1 denotes the state at the lower 
boundary, which in the core collapse case is the NS 
surface or neutrino sphere. The momentum equation, 
eq. (2), in combination with the generic integral equa¬ 
tion, eq. (12), yields 


v s r 2 s (p_u_ - p+v+) = (Pi + pivl)rl - (P + + p+v\)r 2 s 
2 Prdr- j GMpdr. (17) 

In deriving these equations, we have made our first 
major assumption; the solutions between the neutron 
star and the shock are steady. This implies, for example, 
that piV\r 2 = p_v-r 2 , but this does not necessarily 
imply that p_u_ = p+v + . For the sake of brevity, we 
express the right-hand side of Eq. (17) as \F. 



T = (Pi + p x v\)r\ - (P+ + p+v\)r 2 s 

2Pr dr — j GMpdr. ( 18 ) 
Jri 

In the limit of low velocities behind the shock, tF is sim¬ 
ply the overpressure behind the shock compared to hy¬ 
drostatic equilibrium. The utility of using 'F is that it 
incorporates the entire integral solution, including the 
boundary conditions. 

The left-hand side of Eq. (17) must also satisfy the 
momentum jump condition, implying that 

* = (P- + P-vt ~ P+ ~ P+v\)r 2 s . ( 19 ) 

Substitution of this equation into eq. (15) leads to the 
following dimensionless equation for the shock velocity 



v s = -1 + 


\ 


fet) , _ *P 


( 3-1 


(/3 ~ 1 )p+t 


( 20 ) 


s 


where /? is the compression of density across the shock. 
While formally, this expression is correct within our 
given assumptions, its relative complexity hides its util¬ 
ity. With some simple assumptions, we produce a much 
simpler expression between v s and 'F. First, note that 
in the strong shock limit, the ratio of densities, j3, does 
not change as much as the jump in pressure, so that 
the variation of \F dominates the variation in v s . If 
the dynamics of the shock were merely dominated by 
a gamma-law EOS with 7 = 4/3, then the shock com¬ 
pression would be (3 = 7 . However, photodissociation 


of Fe and He nuclei at the shock alters the energetics 
as material passes through the shock. This loss of ther¬ 
mal energy causes a large shock compression. For more 
details on this, see Fernandez & Thompson (2009), for 
example. As a result, the shock compression in most 
core-collapse situations is (3 ~ 9. With such a high shock 
compression, terms such as (3 — 1 ~ /3, which further re¬ 
duces the complexity in eq (20). Therefore, with fairly 
reasonable approximations, eq. (20) reduces to this more 
illuminating expression 

v s ^~i + VTT*, ( 21 ) 

where \F = T /(p+v+r 2 ) is the dimensionless overpres¬ 
sure normalized by the pre-shock ram pressure p+v+r 2 . 
In this form, it is clear that when >F > 0, the shock 
velocity is greater than zero, v s > 0. 

4. VALIDATING THE STEADY-STATE 
ASSUMPTION AND T > 0 WITH 
ONE-DIMENSIONAL PARAMETERIZED 
SIMULATIONS 

Later, we use IF to propose an explosion diagnostic for 
core-collapse simulations, but first we verify in this sec¬ 
tion that \F = 0 during the stalled-shock phase and that 
IF > 0 during explosion. To perform these validating 
tests, we calculate \F in one-dimensional parameterized 
simulations. 

The one-dimensional simulations are calculated us¬ 
ing the new code Cufe, which will be described fully 
in J. W. Murphy & E. Bloor (in preparation). Cufe 
solves the equations of hydrodynamics eqs. (1-3) us¬ 
ing higher-order Godunov techniques. The grid is logi¬ 
cally Cartesian, but employs a generalized metric allow¬ 
ing for a variety of mesh geometries. For this study, 
we use spherical coordinates. The progenitor model 
is the 12 M 0 model of Woosley & Heger (2007), the 
EOS includes effects of dense nucleons around and above 
nuclear densities, nuclear statistical equilibrium, elec¬ 
trons, positrons, and photons (Hempel ct al. 2012). For 
the neutrino heating, cooling, and electron capture we 
use the approximate local descriptions of Janka (2001) 
and Murphy et al. (2013). 

We recall that there are five important parameters 
that describe the steady-state accretion solutions: L u , 
T„, M n g, R u , and M.. Two of these, L„ and T„, we set 
in the parameterized simulations; the rest we calculate 
self-consistently given the equations of hydrodynamics, 
the progenitor structure, and the EOS. The top panel of 
Figure 2 shows the time evolution of Mns, -R„, and 
we highlight the post-bounce phase of a model that does 
not explode. M starts quite high, well above 1 M 0 /s, 


4 We used the SFHo relativistic meanfield variant. 








but then drops down to ~ 0.2 M 0 /s near the end of 
the simulation. At around 500 ms after bounce, A4 has 
a significant drop, which is a result of a density shelf 
advecting through the shock. Throughout this paper, 
you will notice that this significant drop in A4 becomes 
imprinted on many of the explosion diagnostics, includ¬ 
ing the shock radius vs. time in the bottom panel of 
Figure 2. 

When compared to more realistic simulations 
(Melson et al. 2015, see Figure 4), R v evolves very lit¬ 
tle during the time shown in Figure 2. This is a known 
drawback of the standard light-bulb prescription, but 
it does not impact the qualitative conclusions that we 
present later. In more realistic calculations such as 
Melson et al. (2015), the proto-neutron star contracts 
for two reasons. For one, as matter piles onto the PNS, 
the neutron star compresses a little. The primary reason 
that the proto-neutron star contracts, however, is that 
the core neutrino luminosity cools the neutron star. In 
our neutrino model, we omit the diffusive cooling of neu¬ 
tron star. Hence, our core only contracts as a result of 
the added weight of matter. 

However, because each moment may be modeled as 
a successive set of steady-state solutions, this lack of 
contraction does not affect our general conclusion that 
we may calculate an explodability parameter. Later, we 
will calculate an explodability parameter that depends 
upon the five parameters, R v being one of them. Be¬ 
cause the solutions are time independent, the explod¬ 
ability parameter calculation is also time independent. 
Each moment in time has its own explodability value 
that is independent of any other moment. Therefore, it 
does not matter what the evolutionary history of R u is; 
we are able to calculate the explodability parameter for 
any neutrino-sphere radius history. 

To sample a range of explosion timescales, we vary the 
light-bulb neutrino luminosity, L v . The bottom panel 
of Figure 2 shows the resulting shock radii vs. time la¬ 
beled by L v in units of 10 52 erg s -1 . In all models, the 
shock forms at 153 ms after the start of the simulations 
and quickly stalls between ~H50 to ~200 km. The gen¬ 
eral trend is that the higher luminosity models explode 
earlier than lower luminosity models. The lowest lumi¬ 
nosity model does not explode at all. However, it does 
experience a significant outward adjustment and oscil¬ 
lation of the shock as A4 drops significantly around 500 
ms after bounce. The L„ = 2.5 & 3.0 models explode 
during the advection of the density shelf. An obvious 
feature present in the shock radius evolution plot are the 
shock radius oscillations. They are often present in one¬ 
dimensional simulations near explosion (Ohnishi et al. 
2006; Buras et al. 2006b; Murphy & Burrows 2008), and 
Fernandez (2012) suggests that they might be related to 
the advective-acoustic feedback loop that is responsible 




Figure 2. Characteristics of one-dimensional parameterized 
simulations. Top panel: mass-accretion rate (At), mass of 
neutron star (Mns), and neutrino sphere or radius of neu¬ 
tron star (R v ) vs. time for the lowest luminosity simulation. 
Later we compare the integral condition for explosion with 
parameterized simulations. Here, we show the three impor¬ 
tant parameters of the steady-state problem highlighted in 
Figure 1, which are calculated self-consistently in the sim¬ 
ulations. The evolution of these parameters with time is 
similar for all simulations, therefore we only show the lowest 
luminosity simulation here. Bottom panel: shock radius ( R s ) 
vs. time for four one-dimensional parameterized simulations. 
To facilitate comparison with the integral condition, we pa¬ 
rameterize the neutrino heating and cooling. The neutrino 
luminosities are in units of 10 52 erg s -1 . 


for the standing accretion shock instability (SASI) in 
multi-dimensional simulations. 

Our first task is to verify one of the primary results 
of this manuscript, the relationship between v s and the 
integral quantity 4', eq. (21). Figure 3 plots the shock 
velocity normalized by the pre-shock in-fall speed, v s = 
u s /v + , for the L„ = 3.0 model. For comparison, we 
plot the right-hand side of eq. (21), — 1 + \/l + 4/, where 
4/ is calculated directly from the simulation. For the 
most part, the curves agree, validating our derivation 
and assumptions. 

Next, we validate the steady-state assumption. Before 
we do so, we need to consider how time dependence 
might affect eq. (21). In steady-state, the momentum 
integral equation is v s r 2 (p-V- — p+v+) = 'k, but if we 
include the time-dependent term, then the momentum 
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Figure 3. Normalized shock velocity vs. time after bounce. This plot validates the assumptions and derivation of eq. (2 ), 
which relates the shock velocity v s to pressure integral condition 4/. v s is the shock velocity normalized by the the pre-shock 
in-fall speed, v+, and T is overpressure behind the shock in integral form. When 4> « 0, v s ~ 0, and when 'k > 0, v s > 0. This 
validation sets the stage for the explodability condition in Figure . 


integral equation is 

v a rl{p-v- - p+v+) =x + ^, (22) 


where 


X = 



d(pv) 

dt 


dV, 


(23) 


and n is the size of the proto-neutron star, or neutrino- 
sphere radius. Our goal is to see how this time-derivative 
term affects our expression for the shock velocity, there¬ 
fore we use eq. (22) in our derivation for the shock ve¬ 
locity in section 3.2. The equivalent of the final result, 
eq. (21), is 

v fl «-l +Vl + tf+ X, (24) 

where x is just X normalized by p+u^r 2 . 

In some sense, Figure (4) has validated Equation (21) 
in which we ignored the time-derivative term, so that 
it already hinted that x would be small. To be certain 
that x ~ 0, we compare 4» and x in Figure (4). In gen¬ 
eral, (1 + v s ) 2 — 1 = T' + x, but if x ~ 0, then we expect 
(1 + v s ) 2 — 1 ~ 4>. During the steady-state phase, when 
the average shock velocity is zero, we find that both 4> 
and x are small and approximately zero. In fact, their 
small amplitudes nearly cancel. During the initiation of 
the explosion, 'F grows substantially, but x ~ 0- Later, 


when the explosion really takes off, both x and Ik are 
important. Interestingly, we find that assuming steady- 
state is a valid approximation during the steady-state 
phase and the initiation of the explosion, but once explo¬ 
sion commences in earnest, one must clearly consider a 
time-dependent evolution. The focus of this manuscript 
is in deriving a condition for the initiation of explosion, 
however, in which case Figure 3 suggests that eq. (21) 
is a fine place to start. 

In Figure 5, we show 4/ for all of the models that we 
considered. ’F is indeed zero during the steady-state 
phase and 4/ > 0 during explosion. The fact that simu¬ 
lations roughly validate 4f > 0 implies that 4f is a good 
dimensionless explosion diagnostic. 

5. T min : A NEARNESS-TO-EXPLOSION 
CONDITION 

The fact that 4^ > 0 during explosion (Figure 5) sug¬ 
gests that 4> has the potential to be a good explosion 
diagnostic. However, just calculating 4f from the simu¬ 
lations is not the most useful condition; 4* remains near 
zero during the non-exploding phase, and it only devi¬ 
ates from zero while the simulation is exploding. This 
behavior is not a very useful explosion diagnostic, and 
one might as well only use the shock radius. In fact, 
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state equations 



Time after Bounce [s] 

Figure 4. Checking that steady-state is a reasonable as¬ 
sumption in the integral condition. Here we show the over¬ 
pressure, and time-derivative of momentum, Xi as a func¬ 
tion of time after bounce. During the stalled-shock phase, 
the average shock velocity is very nearly zero, and and x 
oscillate about zero as the shock oscillates. Even with these 
oscillations, both and x are quite small. During the initial 
stages of explosion, x ~ 0, and the amplitude of the shock 
velocity is best represented by W alone. Later, as the ex¬ 
plosion proceeds vigorously, both 'I' and x contribute to the 
velocity. Since we are concerned with the revival phase of 
the shock, we find that steady-state is a reasonable approxi¬ 
mation during initiation of the explosion. 


most explosion “conditions” or “diagnostics” to date 
have this problem. They do not provide a nearness-to- 
explosion condition. What we need is a way to translate 
the > 0 condition into a useful nearness-to-explosion 
condition. Fortunately, lends itself to such an explo¬ 
sion diagnostic. 

Our strategy for developing a nearness-to-explosion 
condition is to extract L v , T„, R v , -Mns, and A4 from 
the simulations and calculate the quasi-steady-state so¬ 
lutions for this set of parameters. For a given set of 
parameters, there is a family of solutions to the quasi¬ 
steady-state equations, each with a different shock ra¬ 
dius. Each solution has a value for 'k which may 
be < 0, = 0, or > 0, which corresponds to solutions 
with v s < 0, v s = 0, or v s > 0, respectively. The 
v s = 0 (*F = 0) solution is a quasi-equilibrium solu¬ 
tion. We will show that for each family of solutions, 
there is always a minimum 'k, which we denote v f , m in- 
When tj/min < 0, then v s = 0 solutions exist. How¬ 
ever, when T m i n > 0, the only solutions that exist have 
v s > 0. Therefore, we propose that ’k m ; n is an excellent 
explosion diagnostic providing a nearness-to-explosion 
condition. Figures 6 & 7 illustrate these points. 

Having outlined the strategy, we now provide the 
specifics in calculating ’F m i n . We begin with the steady- 

5 Far from the equilibrium shock radius, these solutions are 
likely incorrect in detail since they correspond to situations that 
are far from a steady-state. However, our approach does not re¬ 
quire accurate models far from equilibrium 


M = Aitr 2 pv , 

(25) 

dv dP (j) 

^ dr + dr ^ r 

(26) 

M de d(r 2 v) 

47rr 2 dr ^ r 2 dr 

(27) 


To solve these equations, we need an EOS to relate 
P in terms of e and p , and we need boundary condi¬ 
tions. For the region that we consider, above the neu¬ 
trino sphere and below the shock, the dominant con¬ 
tributors to the EOS are neutrons, protons, helium nu¬ 
clei (a), positrons, electrons, and photons. Because this 
region is not dominated by dense nuclear physics, the 
EOS is more straightforward and for the most part ana¬ 
lytic. Eventually, we hope to translate the results of this 
manuscript into an analytic solution for the conditions 
for explosion. To facilitate that later goal, we use the 
analytic EOS now. We treat the neutrons, protons, and 
as as an ideal nonrelativistic gas; the positrons, elec¬ 
trons, and photons constitute the relativistic part of the 
plasma, and we consider the positrons and electrons in 
an arbitrary degeneracy. To calculate the relative abun¬ 
dances of neutrons, protons, and as, we assume nuclear 
statistical equilibrium for these three components and 
use the Saha equation to calculate their abundances. 
See appendix B for the details in calculating the EOS 
and the abundances. 

In the spirit of Burrows & Goshy (1993) and others 
(Yamasaki & Yamada 2005, 2006; Yamasaki & Foglizzo 

2008), we solve the steady-state equations, Eqs. (25- 
27), using the Rankine-Hugoniot jump conditions for a 
stalled shock. The jump conditions give the post-shock 
state (subscript —) in terms of the pre-shock state (sub¬ 
script +): 


P-V- 

= P+V +, 

(28) 

V + p _ 

= P+v\ + P+ , 

(29) 

2 

9 


7 V - 

h ~ + Y 

b + h 

+ 

+ 

II 

(30) 


For the pre-shock state, we make the following assump¬ 
tions that enable analytic solutions of the conditions. 
First, we assume that M = A r Kr 2 p + v + is a constant. 
Second, we assume that the star is accreting onto the 
shock at a large fraction of free fall, so the pre-shock 
speed is given by v+ = 2 a v <t>(R s ), where 4>(R S ) is the 
potential at the shock radius, and a v is the fraction of 
free-fall. Third, we assume that the pre-shock Bernoulli 
constant is roughly zero, h+ + v\ — ^>(R S ) ~ 0. Finally, 
we assume a gamma-law relationship for the pre-shock 
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Figure 5. the integral explosion condition, evaluated from the one-dimensional simulations as a function of time. In section 3, 
we derive an integral condition for explosion based upon the momentum equation; we derive that v s > 0 corresponds to 'k > 0. 
To test whether this condition corresponds to explosion, we evaluate for the one-dimensional parameterized simulations 
described in section 4; we find that d/ = 0 during the stalled-shock phase, and > 0 during explosion. These results strongly 
suggest that > 0 is a useful condition for explosion. However, this condition does not by itself predict explosion. It only 
indicates that explosion has occurred. In section 5, we propose a method to turn 4/ > 0 into a useful explosion diagnostic. 


pressure, internal energy, and density: 74 = P/(pe) -f 1 . 
Under these assumptions, the pre-shock pressure is 

P+ = ~ av) 4>{R S ) ■ (3i) 

74-1 

Specifying JVi, a v , 74 , and R s then sets the boundary 
conditions at the shock. When Burrows & Goshy (1993) 
first introduced the critical luminosity, they specified 
one more boundary condition at the base: r = f npdr = 
2/3. This extra condition enables one to self-consistently 
solve for the shock radius that permits a steady-state- 
stalled-shock solution. Yamasaki & Yamada (2005) 
noted that the density at the radius where r = 2/3 is 
almost always the same, so that one can easily replace 
the r condition by specifying a specific density, P 2 / 3 , at 
the neutrino sphere, or neutron star surface. The tradi¬ 
tional way to find the critical curve is to look for solu¬ 
tions that satisfy p\ = P 2 / 3 ', above a curve in luminosity 
and accretion-rate space there are no solutions that sat¬ 
isfy this condition. Burrows & Goshy (1993) interpreted 
this as a critical condition for explosion. However, it has 
never been clear why this condition on the inner density 
should naturally lead to a critical condition for explo¬ 
sion. Now, we show that the solutions above the curve 


indeed only have solutions that correspond to v s >0. 

To do this, we show that a condition on the inner 
boundary, P 2 / 3 /P 1 > 1, is equivalent to Ik > 0 and 
v s > 0. Historically, the method for finding the critical 
curve is as follows. First, for a given R S: one finds “so¬ 
lutions” to the steady-state equations by starting with 
the v s = 0 jump conditions at the shock and integrat¬ 
ing steady-state equations inward to the inner bound¬ 
ary. The resultant partial solution satisfies the govern¬ 
ing equations except that they do not necessarily satisfy 
the inner boundary condition, which means that they 
are not a true solution. One then finds the true solution 
by modifying the shock radius until the partial solution 
also gives the correct inner boundary condition. One 
may use these partial solutions at any shock radius to 
infer what v s would be as a function of R s . 

Understanding the density profile and pressure pro¬ 
files of these partial solutions is the key to connecting 
the old condition on the inner density with a condition 
on v s . Note that ^ in eq. 21 depends upon Pi, p, and p\. 
If one could specify analytic solutions for the pressure 
and density profile, then the condition would be ana¬ 
lytic. In a future paper (J. W. Murphy & J. C. Dolence, 
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in preparation), we will do just that. For now, we use 
the numerical solutions to the steady-state equations to 
propose a semi-analytic solution for the explosion con¬ 
dition. It is our experience that the density profile of 
the partial solution is similar for a given set of param¬ 
eters and R s . In other words, the density profile may 
be written as p{r) = p\z(r). Furthermore, a natural 
dimensionless expression relating the density and pres¬ 
sure is y = P/(p4>) 7 and we find that the shape and 
scale of y mostly depends on the five parameters and 
R s . Therefore, we find that P is also proportional to p\ 
via P = p\y(r)z{r)<j>{r). Substituting these expressions 
for p and P into the expression for 'F (Equation 18), one 
finds that 


^ = Pif ~ P+r 2 s ~ P+v\r 2 s , (32) 

where / incorporates y(r), z(r), and <j>(r) in the first, 
third, and fourth terms of eq. (18). When one inte¬ 
grates the steady-state equations, one naturally finds 
that \F = 0 but the inner density p\ is not necessar¬ 
ily equal to the desired inner boundary, p 2 /3- Because 
p = piz(r), then we may easily substitute p 2 /3 to see 
how 'F changes when we demand that pi = p 2 /3- hr 
terms of pi, the inner density of the partial solution, 
and P 2 / 3 , the desired inner boundary, the dimensionless 
overpressure becomes 


'F 

(P+ + P+v 2 + )r1 


P 2/3 _ 1 
Pi 




(33) 


where the final approximation comes about because the 
pre-shock pressure is typically much lower than the ram 
pressure, P+ -C p+v+- By eq. (21), this approximate 
expression for ’P gives the following expression for the 
shock velocity in terms of the inner densities: 


v s 



(34) 


This presents our final derivation to show that the so¬ 
lutions above the critical luminosity-accretion-rate curve 
have positive shock velocity. Before this derivation, 
Burrows & Goshy (1993) and many others suggested 
that the inability to find solutions that have pi = p 2 /3 
implies explosion. In fact, our derivation implies the 
following explosion condition: 

^ >1—S>>F>0—^v s >0. (35) 

Pi 

To find $ as a function of R s , we use the one¬ 
dimensional simulations to inform the values for A4, 
MnS; Rvi oe „, and 74 , and of course, we set L v and T v to 
the values used for the one-dimensional parameterized 
simulations. Then we find solutions to the steady-state 
equations for a wide range of shock radii. These solu¬ 
tions represent a family of solutions at different R s but 


all with the same parameters. For each solution and as¬ 
sociated R s , we evaluate \F using eq. (33), where we fol¬ 
low Yamasaki & Yamada (2005) and set p 2 / 3 = 7 x 10 10 
g cc -1 . Since v s ss — 1 + y/1 + \F, T naturally shows 
which solutions have v s < 0, v s = 0, or v s >0. The 
v s = 0 solution corresponds to the steady-state-stalled- 
shock solution. 

Figure 6 shows the outcome of this process. Each line 
corresponds to a specific set of parameters and one fam¬ 
ily of solutions. For clarity, we only vary L„ between 
different families of solutions. For low values of L„, all 
three solutions (v B > 0, v s = 0, and v s < 0) are pos¬ 
sible. For low R s , T > 0, implying that v s > 0 and 
the shock would move outward. For larger R s , T < 0 
so that v s < 0 and the shock would tend to move in¬ 
ward. In between, at a very specific shock radius, there 
is a solution that has \F = 0 and v s = 0. This solu¬ 
tion represents an equilibrium solution. While solutions 
far from equilibrium are not steady and therefore are 
poorly represented by the steady-state equations, our 
approach does not rely on the quantitative accuracy of 
this representation. For a visual aide to these concepts, 
see Figure 7. 

Note that for high values of there are no solutions 
for which 'F min < 0. In these cases, all solutions corre¬ 
spond to « 5 >0. There is a very specific set of param¬ 
eters for which 'P m in = 0. The locus of such parame¬ 
ters defines a critical hypersurface that generalizes but 
encompasses the critical neutrino luminosity condition 
of Burrows & Goshy (1993) (see Section 6 for more on 
this). 

To validate that 'F m ; n satisfies the desired qualities 
of an explosion diagnostic, we plot ^mm in Figure 8 for 
the one-dimensional parameterized simulations. The pa¬ 
rameters that we set by hand in each simulation are L„ 
and T v ; we keep T v set at 4 MeV, and vary L At each 
time, we extract the other important parameters Mns, 
I?„, and A4, which are calculated self-consistently in the 
simulation. We then calculate the family of steady-state 
solutions for that set of parameters and find the value of 
'I'min- Figure 8 shows that \F m i n has the desired qualities 
of an explosion condition. Before explosion, ’Fmin < 0 
and during explosion, Tmin > 0. Moreover, we suggest 
that the value of 'Fmin before explosion provides a useful 
metric for how far away the simulation is from explosion. 

6. COMPARISON WITH OTHER EXPLOSION 
CONDITIONS 

In this section, we briefly compare the integral 
explosion diagnostic, 'Fmin with three other con¬ 
ditions: the critical neutrino luminosity condition 

(Burrows & Goshy 1993), a timescale ratio condition 
(Janka & Keil 1998; Thompson 2000; Thompson et al. 
2005; Buras et al. 2006a; Murphy & Burrows 2008), and 
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Figure 6. The family of steady-state solutions represented by their values of \l/( R a /R „) for three different neutrino luminosities 
but all else being equal. For some parameters, W may be negative, zero, or positive. The solution with T = 0 (u 3 = 0) 
corresponds to the equilibrium solution. For each family, there is a solution with a minimum 'll, which we denote Tmin. When 
min < 0, a stalled-shock ( v s = 0) solution exists, but when 'I'mm > 0, only solutions with v s > 0 exist. We propose that 
'll min > 0 is an excellent explosion diagnostic (see Figure ). 


the antesonic condition (Pejcha & Thompson 2012). In 
this manuscript, we mostly compare the effectiveness of 
these conditions with the integral condition; in a forth¬ 
coming paper (J. W. Murphy & J. C. Dolence, in prepa¬ 
ration), we show that one can actually derive all three 
conditions from the integral condition by making various 
further approximations. For now, we simply compare 
the conditions. 

Of these three, the most closely related condition is 
the critical luminosity condition of Burrows & Goshy 
(1993). In fact, a primary motivation in deriving the 
integral condition is to derive a condition that shows 
that the solutions above the critical neutrino luminosity 
curve correspond to v s > 0. In section 3, we derived the 
integral condition for v s > 0, and now, in the bottom left 
panel of Figure 9, we show that this same integral con¬ 
dition reproduces the critical neutrino luminosity curve 
of Burrows & Goshy (1993). 

To reproduce the critical neutrino luminosity curve 
in the lower left panel of Figure 9, we first fix three 
of the five important parameters of the problem, 1?„, 
Mns, and T„. In other words, we restrict the dimen¬ 
sionality of the integral condition to the L V -J14 plane. 
Then we find the solutions to the steady-sate equations 


(eqs. 25-27) and select the solutions in the L u -A4 plane 
that have d'mm = 0. The locus of these specific solu¬ 
tions forms the solid curve that we show in Figure 9. 
Above this curve, 'Pmm > 0 and therefore all solutions 
have v s > 0. Below this curve, 'P min < 0, and there¬ 
fore a steady-state stalled-shock solution exists. This 
curve is exactly the same critical neutrino luminosity 
curve that Burrows & Goshy (1993) and others have 
derived. The difference is in the way that it is de¬ 
rived. Burrows & Goshy (1993) solved the steady-state 
equations looking for solutions for which the shock and 
r = 2/3 conditions are satisfied. There is no statement 
about the behavior of v s above the curve. Rather than 
highlighting r, we instead use the ’F m i n = 0 condition, 
that one readily sees that v s > 0 above the curve. 

In addition to showing that v s > 0 above the critical 
neutrino luminosity curve, we also show that ’F m i n > 0 
is a more general explosion condition than the criti¬ 
cal L„-A4 curve. In fact, in the important parame¬ 
ters, 'F m ; n = 0 represents a critical hypersurface, and 
the critical L v -Jvl curve is just one slice of this more 
general critical condition. To help illustrate this, we 
use in Figure 9 the ’F m in = 0 to derive other critical 
curves, where each is just a different slice of the criti- 
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Figure 7. This is an illustration to demonstrate the two 
regimes seen in Figure 6. If ^min < 0, then there are so¬ 
lutions with v s < 0, v a = 0, and v 3 > 0. In general, the 
solution with v s = 0 is an equilibrium solution, and as long 
as it exists, simulations (and presumably Nature) settle on 
this steady-state stalled-shock solution. When Tmin > 0, the 
only solutions are those in which v s > 0. Therefore, T m i n 
seems to be a natural condition for explosion. 

cal hypersurface. Each critical curve is constructed in 
the same way as the L u -A4 critical curve, except that a 
different set of three parameters is fixed. In all panels, 
T v is fixed to 4 MeV; each panel shows the values for 
the other two parameters. In the past, the literature 
has focused on the L„- A4 critical curve, mostly because 
that is how this useful critical curve was discovered and 
presented (Burrows & Goshy 1993). However, Figure 9 
clearly shows that there are many other critical curves, 
one for each pair of the five parameters. 

By focusing on one critical curve, one runs the danger 
of overemphasizing two of the parameters and underem¬ 
phasizing the other three parameters. For example, the 
L u -M critical curve suggests that one should pay atten¬ 
tion to the evolution of L u and A4. When it comes to 
the equations, however, there is nothing special about 
L„ or Al; all parameters are potentially important. De¬ 
pending upon the evolution of the parameters, it might 
therefore be more instructive to consider the evolution 
of, for example, R v and A4. 

The Tmin = 0 critical condition suggests another way 
to view the nearness-to-explosion that does not overem¬ 
phasize any one parameter. The five parameters define 
a five-dimensional hyperspace, and the 4' m i n = 0 condi¬ 


tion defines a hypersurface within the five-dimensional 
space. Below the 4'min = 0 hypersurface, steady-state 
solutions exist, and above this hypersurface we show 
that the solutions are likely explosive. As a simulation 
evolves, it will trace out a path in this 5-dimensional 
hyperspace, and in general, it will move toward or away 
from the 4 , m ; n = 0 hypersurface. Therefore, a general¬ 
ized distance to this surface would be an excellent mea¬ 
sure of nearness-to-explosion. In a future paper, we will 
present discussions of this generalized distance. For now, 
we merely introduce this new concept, and propose that 
since the value of 4'min defines the hypersurface, it is a 
good proxy for the generalized distance to the bounding 
hypersurface. This definition is useful for two reasons. 
One, since the hypersurface is defined by 4'min = 0, we 
are able to define the critical condition by one dimen¬ 
sionless parameter and not five. Two, this dimensionless 
parameter is intimately related to a fundamental ques¬ 
tion of CCSN theory. 4' is proportional to v s , and a 
fundamental question of CCSN theory is “what are the 
conditions for v s > 0?” Hence a single dimensionless 
parameter, 4'min, that defines the condition for explo¬ 
sion is also intimately related to an important question 
in CCSN theory. 

To help illustrate the critical hypersurface, Figure 10 
shows two slices of this hypersurface. Since it is difficult 
to visualize a five-dimensional hypersurface, we instead 
show two slices of the hypersurface. In the top panel of 
Figure 10, we fix T v to 4 MeV and R u = 63 km. The 
resulting surface is a two-dimensional surface defined by 
three dimensions. In each of the axis planes, we also 
show critical curves, with each curve corresponding to a 
specific value of the third axis. The bottom panel shows 
an equivalent slice, but this time the slices are at T v = 4 
MeV and = 3 x 10 52 erg/s. While it is impossible 
to show the full hypersurface, we hope to illustrate with 
these two-dimensional slices that the critical condition 
is not only a critical curve, but a critical hypersurface 
defined by 4'min = 0. Below the hypersurface, 4'min < 0, 
and steady-state solutions exist; above the hypersurface, 
4’min > 0 and all solutions likely have v s > 0. 

Another common class of explosion metrics are the 
timescale ratios in which one compares an advection 
timescale (r a d v ) to a heating timescale (rh ea t)- Roughly, 
neutrino heating might be expected to significantly mod¬ 
ify the internal energy of the accreting material only 
when Th ea t ;$ T adv- Based on this, many have used 
the ratio t^/t h ea t as a diagnostic to indicate how 
close a model is to explosion for a given snapshot 
(Janka & Keil 1998; Thompson 2000; Thompson et al. 
2005; Buras et al. 2006a; Murphy & Burrows 2008). 
While this may seem sensible at face value, the timescale 
ratio arguments all neglect important aspects of the 
CCSN problem — cooling that occurs below the gain 
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Figure 8. Nearness-to-explosion condition for one-dimensional hydrodynamic simulations: 'F m m vs. time. For each time, 
we extract the important parameters of the problem (L v , T v , R„, Mas, and M) from the one-dimensional hydrodynamic 
simulations and calculate the family of steady-state solutions. 'I'min represents the minimum possible value of 'P over that 
family (see Figure 6). While 'I'min < 0, there exists a stalled accretion shock solution (4> = 0 and v s =0). When 'I'min > 0, 
all possible solutions have v 3 > 0, which we associate with explosion. Empirically, we find that before explosion 'I'min < 0 and 
during explosion ’I'min > 0. Most importantly, the value of 'I' m ; n offers an excellent indication of nearness-to-explosion. 


region, for example. Meanwhile, urphy & Burrows 
(2008) and others more recently (Dolence et al. 2013; 
Takiwaki et al. 2014; Hanke et al. 2012) have shown 
that the turbulent dynamics of multidimensional lead 
to longer advection times on average, speculating that 
this may be an important effect vis-a-vis ID vs. multi-D. 

There are many approximate definitions of r ac j v and 
Theat (Janka & Keil 1998; Murphy & Burrows 2008; 
Fernandez 2012) . Early definitions of the advection 
time include r a d v = f T " hock dr/v r , but the most recent 

J ' gain ' 

definitions use the mass in the gain region and the mass- 


accretion rate 


Fulv 


M. 


gam 


M 


(36) 


For the rest of this paper, we adopt this most recent def¬ 
inition. For the heating timescale, the generic approach 
is to compare a total energy with an integrated heating 
rate Th ea t = E/Q. Generically, many define the heating 
rate as Q = J gain (H — C)pdV. The total energy has sev¬ 
eral definitions, and since none are derived from a firm 
explosion condition, all are arbitrary. We highlight and 
use one condition; we consider the total internal energy 
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Figure 9. “Critical curves” defined by the 'h m i I1 = 0 condition. These plots show that the 'h m i n = 0 is another way to define 
the critical curve, the solutions above the critical curve have positive shock velocity, and the L v -M curve is merely a slice of 
a critical liypersurface. Each panel represents a two-dimensional slice through the five-dimensional parameter space. In all 
panels, T„ is fixed to 4 MeV. The two other fixed parameters are labeled in each panel. The color map shows the value of 'k m i n 
as a function of the two remaining parameters, and the solid black line shows where TnUn = 0. In the lower left panel, we 
show that the ^min = 0 condition produces the critical L„ — M curve of (19 ). Note that above this curve 

\kmin > 0, which implies that solutions above the curve have v s > 0. In the context of our assumptions, this suggests that the 
solutions above the critical curve are explosive, as un ( ) initially suggested. Finally, with these panels, we 

show that the critical neutrino luminosity curve is a slice of a more general critical hypersurface. Most importantly, this critical 
hypersurface is represented by one dimensionless condition, \k m i n = 0. 


in the gain region E = J gain £p dV. The choice of both 
the energy scale and the region are completely arbitrary 
and not derived from an actual explosion condition. In 
a forthcoming paper (J. W. Murphy & J. C. Dolcnce, in 
preparation), we will derive a timescale ratio condition 
from the integral condition, but for now we just compare 
to the arbitrary definitions from the past. We pick one 
example, the total internal energy in the gain region, 


Figure 11 compares our derived integral condition, 
Ikmin, with the approximate timescale ratio condition, 
T a dv/Theat- As one might expect, the rough heuristic 
condition does not perform as well at predicting explo¬ 
sions. Indeed, r a d v /Theat is 0(1) before explosion and in¬ 
creases dramatically after explosion. However, the ratio 
has serious problems as an explosion metric. For exam¬ 
ple, one might be tempted to conclude that the sharp 
upward trend in this ratio would be a good indicator 
of explosion. However, careful inspection of T a d v /Theat 


and the shock radii in Figure 11 shows that T a dv/iheat 
just follows R s . Therefore, using r a dv/iheat is not much 
more useful than R s . One might be tempted to use the 
simulations to calibrate a “critical” value for 'r a d v /'Tieat- 
However, the results in the literature as well as those 
shown here indicate that no such “critical” value exists 
(Muller et al. 2012b; Hanke et al. 2012; Dolence et al. 
2013). 

The final comparison is between 'F m in and 
the antesonic condition of Pejcha & Thompson 
(2012); see Figure 12. Similar to our motivation, 
Pejcha & Thompson (2012) were motivated to explain 
the origin of the critical luminosity curve and possibly 
derive a more general condition. They took a slightly 
different tack. Pejcha & Thompson (2012) considered 
the family of accretion and wind solutions and hypoth¬ 
esized that explosion occurs at the intersection between 
steady-state accretion solutions and steady-state wind 
solutions. By considering isothermal profiles, they were 
able to define an analytic condition that divides these 
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Figure 10. Critical hypersurface for explosion. The 41 m i n = 
0 condition defines a hypersurface in the five-dimensional 
parameter space. Below the hypersurface, ’I'min < 0, and 
steady-state solutions exist; above the hypersurface, 4' m i n > 
0. Therefore, only solutions with v s > 0 exist. It is difficult 
to visualize the full hypersurface; therefore we present two 
slices that are themselves two-dimensional surfaces. In both 
plots, we fixed T- to 4 MeV. In the top plot, we also fixed 
R v to 63 km, and in the bottom plot we fixed L„ to 3 x 10 52 
erg s' 1 . We also show “critical curves,” for which we specify 
the value for the third dimension. 

solutions. For c^/v g SC > 3/16, only wind solutions 
are possible, which they attributed to the transition 
between steady-state accretion and explosion. Given 
the great deal of neutrino heating and neutrino cooling, 
the post-shock profile is certainly not isothermal. 
Therefore, they were unable to derive a similar analytic 
condition in a more realistic CCSN context. Instead, 
they solved the steady-state equations and empirically 
searched for a similar condition that would be appropri¬ 
ate for the core-collapse case. The resulting explosion 
condition that they proposed is the antesonic condition, 
max(c s /r eS c) 2 > 0.19, where max() denotes the max¬ 
imum in between the NS and the shock. Empirically, 
our one-dimensional simulations indicate that a value of 
~0.19 is about right. In multidimensional simulations, 
Muller et al. (2012a) and Burrows (2013) found a much 
higher value of ^0.3. However, given that the antesonic 
condition was derived and discussed in the context of 
one-dimensional spherical symmetry, it is no surprise 
that the multidimensional simulations do not match. 
After all, we derive an integral condition in the same 



Figure 11. A comparison of the explosion diagnostic, 
41 min > 0, and a timescale ratio condition, Tadv/rheat- The 
top plot shows 4/ m i n and the bottom plot shows the ratio of 
the advection time (r a dv) to the heating time (rheat). The 
advection timescale is defined as r a dv = M ga i n /.M, where 
Mg a in is the mass in between the gain radius and the shock. 
We adopt a heating timescale that is a ratio of energy and 
heating rate in the gain region: Th ea t = E/Q gain, where 
Qgain = (H — C)pdV, and E = J gliin spdV is the total 

internal energy in the gain region. The r a d v /'rheat is an ap¬ 
proximate heuristic explosion diagnostic. The choice of the 
terms to dehne the timescales and the region over which to 
integrate are arbitrary. The timescale ratio is 0( 1) but it 
shows no sign of a reliable “critical” value or utility in pre¬ 
dicting explosion. In fact, it is only as good as using the 
shock radius as an explosion diagnostic. On the other hand, 
41 min works quite well as an explosion diagnostic. 

spherically symmetric context. In a forthcoming paper, 
we rederive the integral condition including turbulence, 
and find that the new integral condition is consistent 
with simulations. We suspect that one could do the 
same for the antesonic condition. However, to use the 
antesonic condition in practice, one must measure the 
“critical” condition in simulations ex post facto. While 
the antesonic condition, which is analytically derived 
in the isothermal case and empirically derived in the 
nonisothermal case, is imperfect, Figure L2 shows that 
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it fairs a little better than the timescale ratio condition. 

While the antesonic condition fares slightly better 
than the timescale ratio at predicting explosion, it still 
remains less favorable compared to the integral con¬ 
dition (middle panel of Figure 12). For one, like the 
timescale ratio, the dramatic upward trend in the an¬ 
tesonic ratio that might be attributed to explosion ac¬ 
tually occurs well after explosion; this means that one 
can rule out using the upward trend as an explosion 
diagnostic. One of the useful features of the integral 
condition, *F m in > 0 is that one can tell even before ex¬ 
plosion that a model is unlikely to explode. There are 
two features that enable this: 1) the closeness of 'Fmin 
to zero, and 2) the trend toward or away from zero; for 
example, 'Fmin for the = 3.0 model inarches toward 
explosion, while 'Fmin for the L v = 2.0 model marches 
away from explosion. Similar to point 1) for the integral 
condition, the antesonic condition is indeed closer to the 
“critical” value before explosion. However, even the an¬ 
tesonic diagnostic for the L„ = 2.0 marches toward the 
“critical” value even though it never explodes. Further¬ 
more, the antesonic diagnostic overshoots the “critical” 
value by ~10% of the “critical” value. This may seem 
like a triumph of estimating the “critical” value, but 
the dynamic range of the antesonic diagnostic is only 
~20% of the “critical” value, making the estimate only 
~50% accurate when considering the dynamic range of 
the antesonic diagnostic. 

7. DISCUSSION AND CONCLUSIONS 

For many years, the critical neutrino luminosity curve 
has been the most impactful model for understanding 
successful explosion conditions. However, no one had 
demonstrated analytically that the solutions above the 
curve are explosive, nor had anyone successfully used 
the critical neutrino luminosity as an explosion diagnos¬ 
tic for core-collapse simulations. In this manuscript, we 
took steps to show that the solutions above the criti¬ 
cal neutrino luminosity have v s > 0, and we used this 
derivation to propose an explosion diagnostic for core¬ 
collapse simulations. When we began this project, our 
primary goal was to merely show that the solutions 
above the critical neutrino luminosity curve are explo¬ 
sive. We have partially done so by deriving an integral 
condition, 'L m i n > 0, for v s > 0 and showing that the so¬ 
lutions above the critical neutrino luminosity curve have 
v s > 0. Along the way to deriving this integral condi¬ 
tion, we discovered that the neutrino luminosity curve 
is a projection of a more general critical hypersurface. 
Most strikingly, this critical hypersurface is represented 
by a single dimensionless parameter, ’Fmin > 0, and it 
promises to be a useful explosion diagnostic for core¬ 
collapse simulations. 

In deriving the integral condition <F > 0, we made two 


simple but profound approximations: (1) we considered 
steady-state solutions, and (2) we considered that v s > 0 
corresponds to explosion. With these two approxima¬ 
tions, we derived an integral condition for v s > 0 and 
recast it in a dimensionless form, *F > 0. We verified 
this integral condition using one-dimensional parameter¬ 
ized simulations and found that during the stalled-shock 
phase 'F ft: 0 and during explosion > 0. Although the 
comparisons with simulations successfully verified the 
integral condition, they also demonstrated that simply 
calculating *F from the simulations is not a useful explo¬ 
sion diagnostic. Because the simulations always find a 
stalled-shock solution before explosion, 'F is always zero 
before explosion. Such a diagnostic fails to indicate a 
distance to explosion. If, on the other hand, we extract 
the important parameters from the simulations and cal¬ 
culate all of the possible steady-state solutions, then a 
better diagnostic emerges. 

The steady-state solutions fall into one of two broad 
categories: in one category, a stalled-shock solution ex¬ 
ists, and in the other, only v s > 0 solutions exist. To un¬ 
derstand the origin of these categories, one must first un¬ 
derstand that there are five important parameters that 
characterize each steady-state solution: L v (neutrino 
luminosity), T„ (neutrino temperature), R u (neutrino- 
sphere radius or proto-neutron star radius), Mns (proto¬ 
neutron star mass), and M (accretion rate). For a fixed 
set of parameters, there is a family of solutions, and 
each solution in this family has a different value of *F. 
In general, 'F may be negative, zero, or positive, and be¬ 
cause *F oc v s , each solution may have v s < 0, v s = 0, or 
v s > 0. If a v s = 0 solutions exists, it represents the pre¬ 
ferred quasi-equilibrium solution. Generically, for each 
family of solutions, there is a minimum 'F that we denote 
^F m i n , and it is this parameter that determines whether 
a solution belongs to the stalled-shock category or to the 
explosive category. If 'Fmm < 0, then the family of solu¬ 
tions has solutions with v s < 0, v s = 0, and v s > 0 with 
the v s = 0 being the preferred quasi-equilibrium solu¬ 
tion. If on the other hand, 'Fmin > 0, then only v s > 0 
solutions exist. If we attribute v s > 0 to explosions, 
then \F m i n is a useful explosion diagnostic. 

We calculated 'Fmin for several one-dimensional pa¬ 
rameterized simulations (see Figure 8) and find that the 
time evolution of ’F m in makes an excellent explosion di¬ 
agnostic. From the time that we track 'Fmin, one can 
tell whether a simulation is near explosion. For those 
that are near explosion, ’F m in tends to march toward ex¬ 
plosion throughout much of the simulation. The model 
that does not explode shows a general evolution of 'F m in 
away from explosion. Next, we need to explore 'Fmin in 
the context of multidimensional and self-consistent neu¬ 
trino radiation hydrodynamic simulations to see if these 
promising behaviors remain. 
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In section 6 and Figures & 12, we compared this 
new integral condition with explosion conditions that 
have been defined previously. Specifically, we compared 
’J'min > 0 to a timescale ratio condition and the an- 
tesonic condition. In contrast to the integral condition, 
the Tadv/^heat condition is heuristic, only good to order 
unity, shows no obvious “critical” value, and is no better 
than simply using R s as an explosion diagnostic. The 
antesonic condition fares better than the timescale ratio 
condition. Even so, it lacks the accuracy and predictive 
attributes of the integral condition. 

To derive the integral condition, we made two approx¬ 
imations that must be investigated further: we assumed 
that v s > 0 corresponds to explosion, and we assumed 
steady-state throughout our derivation. These are two 
very strong assumptions that could have invalidated our 
derivation. For example, v s > 0 does not necessarily 
equate to explosion, especially when we relax the steady- 
state condition. Oscillations of the shock are one exam¬ 
ple in which v s > 0 may manifest but not lead to explo¬ 
sion. In order to make progress, we decided to see where 
these approximations led. In the end, we were able to 
derive an integral explosion diagnostic that successfully 
describes the explosion conditions of parameterized one¬ 
dimensional simulations. However, we did not show that 
once 'Fmin > 0 and the shock begins to move outward 
that it continues to move outward. In some sense, one 
expects the shock to continue to move out, because once 
'J'min > 0 the only steady-state solutions are solutions 
with v s > 0, but this does not preclude the possibility 
of a dynamic (non-steady-state) solution that is oscilla¬ 
tory. A next step in confirming our derivation is to show 
that once ’F m in > 0, the dynamic solution is a solution 
in which v s > 0 continues. 

We found the integral condition to be a success¬ 
ful explosion diagnostic for one-dimensional parameter¬ 
ized explosions, but we suspect that it can do much 
more. We suspect that the integral condition may 
be an excellent explosion diagnostic for self-consistent 
three-dimensional CCSN simulations. One may even be 
able to predict whether a certain progenitor model will 
explode without even performing core-collapse simula¬ 
tions. This hope will only be realized with further model 
developments, replacing parameters currently measured 
from simulations with parameters calculated by other 
means. 

Before we can even pursue such bold endeavors, how¬ 
ever, we must adapt the integral condition to include 
the appropriate physics. For one, we derived the in¬ 
tegral condition using Newtonian gravity; general rela¬ 
tivistic considerations are important, therefore we will 
need to take the straightforward steps in deriving the 
condition in GR. Second, we need to incorporate a more 
self-consistent neutrino heating and cooling. Third, we 


need to incorporate multidimensional effects. Turbu¬ 
lence seems to reduce the critical neutrino luminosity 
for explosion; we suspect that one can easily use a tur¬ 
bulence model to derive a new integral condition for ex¬ 
plosion including turbulence. 

In summary, we derived an integral condition for ex¬ 
plosion and verified it with one-dimensional parameter¬ 
ized simulations. When combined with simple steady- 
state models, it suggests a new explosion diagnostic, 
lk m i n , that we argue may be a more predictive mea¬ 
sure of a models explodability than other diagnostics. 
Finally, we point out that our integral formulation can 
be extended with better physics, and we are hopeful 
that this approach may prove useful in disentangling the 
complicated interactions of various physical effects and 
in understanding the mechanism of CCSNe in Nature. 
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APPENDIX 

A. STEADY-STATE EQUATIONS 

Identifying the steady-state solution with the minimum dt is key in developing the explosion diagnostic in section 5 
and Figure 8. In this appendix, we therefore present the steady-state equations that we used to obtain the family of 
solutions. 

To accommodate the analytic EOS described in the next appendix, we recast the steady-state equations in terms of 
the natural independent variables. Because the analytic EOS in appendix is most naturally written as a function of 
density and temperature, P = P{p,T ), we rewrite the steady-state equations as differential equations for v, p , and T. 
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where P p = (dlnP/dlnp)T is the partial derivative of P with respect to p at constant T. Equivalently, Pt = 
(91nP/91nT) p , e p = (91ne/91n p)r, and er = (9 In e/9 In T) p . 


B. EQUATION OF STATE 

For the one-dimensional parameterized simulations, we use the tabulated EOS provided by Hempel et al. (2012). 
The microphysics includes a distribution of nuclei that satisfy nuclear statistical equilibrium; the individual com¬ 
ponents are, broadly, a dense nuclear component, ideal gas for the nucleons and isotopes, photon gas, and rel¬ 
ativistic electrons and positrons with arbitrary degeneracy. The tabulated EOS and driver are available at 

http://www.stellarcollapse.org/equationofstate. 

For the steady-state solutions we use an analytic EOS that does remarkably well in reproducing the microphysics 
in the region between the neutrino sphere and the shock. Bethe (1990), Janka (2001), and Fernandez & Thompson 
(2009) were useful guides in developing this analytic EOS. First, we assume nuclear statistical equilibrium for three 
species only: neutrons, protons, and as; their respective mass fractions are X n , X p , and X a . Conservation of baryonic 
mass implies 

X n + X p + X a = 1, (B4) 


and conservation of charge implies 

+ ix a = Y e (B5) 

where Y e is the number of electrons per baryon. In our steady-state solutions we simply set Y e = 0.5. In NSE, the 
Saha equation provides the remaining equation to find a solution for the abundance of these three species. 
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and Q a = 28 MeV is the binding energy of the a. 

We construct the pressure and internal energy under the assumptions that the photons, positrons, and electrons are 
relativistic and the partial pressure due to the nucleons and as is given by the ideal gas law. In addition, we consider 
the electrons and positrons with an arbitrary degeneracy 77 = p e /(kBT). The expression for the degeneracy parameter 
is 
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Figure 12. The 'h m i n explosion diagnostic in compari¬ 
son with the antesonic condition of 
(2012). For an isothermal flow in one dimension, 
Pejcha & Thompson (2012) show that there are no steady- 
state accretion solutions that match the pre-shock flow for 
Ct/Ubsc > 3/16. By analogy, they propose that a simi¬ 
lar condition exists for the nonisothermal situation of the 
core-collapse problem. By numerically exploring the steady- 
state solutions in one-dimensional approximations, they pro¬ 
pose that max(Ca/Ue SC ) > 0.19 divides steady-state ac¬ 
cretion from explosion. In multidimensional simulations, 
Muller et al. (2012a) and Dolence et al. (2013) empirically 
find that max(c^/«e SC ) > 0.3 is a better “critical” value. The 
bottom panel of this figure shows that the one-dimensional 
“critical” value is a decent indicator of explosion. It certainly 
performs better than the timescale ratio seen in Figure 
Still, the 'F m in > 0 condition has some practical advantages. 
For one, $ m i n > 0 is an integral condition derived from the 
supposition that v 3 > 0 corresponds to explosion. Second, 
the dynamic range of \E , m i n illustrates more clearly which 
simulations are far from explosion. Third, while the low¬ 
est luminosity model appears to always approach explosion 
for the max(Cs/i>es C ) condition, the ^min condition clearly 
shows that the L v = 2.0 model moves away from explo¬ 
sion, except when the density shelf accretes through the 
shock, but then again it moves away from explosion. Fourth, 
max(Cg/Ue SC ) > 0.19 does not predict explosion: the actual 
antesonic condition overshoots the proposed critical value by 
10 % before it actually explodes, which is half of the full dy¬ 
namic range of the antesonic condition (~ 20%); when the 
antesonic condition looks as if it indicates explosion, it does 
so well after the initiation of explosion. 
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The pressure and internal energy may both be divided into partial pressures due to relativistic ( R) particles and 
nonrelativistic (NR) particles 

P = Pr + Pnr , (B9) 

and 

£ = £r+£nr- (BIO) 


The partial pressure due to the relativistic constituents is 

1 ( k B T ) 4 / 11tt 2 2 rj 4 

R 12 (He) 3 V 15 + V + 7T 2 ) ’ 
and the partial pressure due to the nonrelativistic constituents is 


The resulting internal energy is 
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where the last term is there to account for the transfer of binding energy per nucleon from X a to the gas. 
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